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A rational approach to identify polyatomic molecules that appear to be promising candidates for 
direct Doppler cooling with lasers is outlined. First-principle calculations for equilibrium structures 
and Franck-Condon factors of selected representatives with different point-group symmetries (in¬ 
cluding the chiral non-symmetric group Ci) have been performed and high potential for laser-cooling 
of these molecules is indicated. 
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Introduction Cold molecules are currently heavily 
studied and serve as particularly interesting targets for 
spectroscopy mn]. A wealth of applications of cold 
molecular systems for high-precision spectroscopy includ¬ 
ing the search for violation of fundamental symmetries 
[311], dark matter |S] and variation of fundamental con¬ 
stants |6| as well as controlled cold chemistry |7], quan¬ 
tum computations (see e.g. review |8|) and many oth¬ 
ers, have been considered recently. A special issue of 
Chemical Reviews [5] contains comprehensive reviews of 
some popular techniques for obtaining samples of cold 
and ultra-cold molecules, among which are photoassoci¬ 
ation from cold atoms, Stark ans Zeeman deceleration 
of the molecular beams and buffer-gas cooling. Presently 
the stimulated Raman adiabatic passage (STIRAP) tech¬ 
nique used to transfer Feshbach molecules to the rovi- 
bronic ground state provides lowest temperatures and 
highest densities of diatomic molecular gases (see e.g. |9] 
and recent developments in [Mill]). 

Considerable progress has been reported for direct 
Doppler cooling of diatomic molecules with lasers nans]- 
The latter method looks especially promising, as it is 
intimately connected with well-developed techniques for 
cooling and trapping of atoms. Thus, in principle, diverse 
advances in atomic Doppler cooling can be transferred 
to the related schemes for molecules. For example the 
magneto-optical trapping widely used with cold atoms 
was recently demonstrated for YO [13] and SrF M- 

As of yet, however, only diatomic molecules could 
successfully be Doppler cooled with lasers as molecules 
amenable for Doppler cooling must feature Franck- 
Condon (FC) factors close to unity for transition be¬ 
tween those vibronic states involved in the closed transi¬ 
tion cycles (usually the ground and first excited electronic 
state). This is connected with the fact that Doppler cool¬ 
ing requires absorption and reemission of thousands to 
millions of photons per atom or molecule involved in cool¬ 
ing. The reemission process creates a dissipative force 
for species counterpropogating the laser beam, whereas 
working levels for those species copropogating the laser 
beam are shifted out of resonance because of the Doppler 


effect. Thus, closure of the cooling transition is crucial 
for Doppler cooling. In general, such closed transitions do 
not exist in molecules due to non-selective spontaneous 
decay from the selected rovibronic level of the excited 
electronic state to a multitude of rovibrational levels of 
the electronic ground state. The probability of such a 
process is proportional to the FC factor between the cor¬ 
responding vibrational levels. As it was first displayed 
by DiRosa in 2004 [TS], there exist diatomic molecules 
with the sum of a few largest FC factors (for transitions 
originating from a given vibrational level of the excited 
electronic state) being very close to unity. Peculiari¬ 
ties of the electronic structure in systems possessing such 
quasi-diagonal FC matrices have been identified in Ref. 
13 Almost parallel Born-Oppenheimer potential energy 
curves of the electronic ground and excited state (which 
results in diagonality of the FC matrix) are expected in 
diatomic molecules with one valence electron over closed 
shells, when essentially this electron undergoes transi¬ 
tion between (mainly) non-bonding orbitals upon elec¬ 
tronic excitation. As electronic transition between non¬ 
bonding orbitals influences the chemical bonds (formed 
by electrons on bonding orbitals) typically only weakly, 
the internuclear potential stays practically unchanged for 
ground and excited electronic states. The influence of 
the non-bonding electron is even smaller, when due to 
e.g destructive s-p hybridization the electron is shifted 
away from the bonding region (see next section). 

We pointed out [16] that also in polyatomic molecules 
with a valence electron undergoing transitions between 
non-bonding orbitals quasi-closed transition loops should 
be present, which opens up an avenue to cold sam¬ 
ples of larger molecules. In this Letter we provide ex¬ 
amples of polyatomic molecules that may be expected 
to have quasi-diagonal FC matrices for vibronic tran¬ 
sitions between electronic states and are thus poten¬ 
tially amenable for direct cooling with lasers. Classes of 
molecules and explicit examples together with possible 
routes for generating proper molecular structures were 
originally reported on the Conference for Cold and Con¬ 
trolled Molecules and Ions (Monte Verita (Switzerland), 
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2014). 

General consideration A straight-forward way to 
identify polyatomic molecules with electronic properties 
similar to those of diatomic molecules is to apply sub¬ 
stitution: e.g. in the series of MF open-shell diatomic 
molecules (class I molecules, see [3] and below in this 
section), where M is a metal atom from group 2 of the 
Periodic Table (Be-Ra), one can substitute fluorine for 
a pseudohalogen such as CN, NC, SCN or functional 
groups such as OH or CH 3 . The pseudohalogens are 
polyatomic functional groups whose chemical properties 
resembles those of halogens. Particularly the bonding sit¬ 
uation in molecules containing pseudohalogenic groups is 
close to their halogen analogs, a fact known in chemistry 
for almost a century (see e.g. m)- From the electronic 
structure point of view one might expect that the oc¬ 
cupation pattern of the leading electronic configurations 
for the ground and energetically lowest excited state of 
the resulting compounds differs only for the unpaired 
electron, which occupies then different (essentially) non¬ 
bonding orbitals, similar to the original MF molecule (see 
also [IH] for a more detailed discussion on the expected 
structure for calcium derivatives). For some of the above- 
mentioned compounds we also benefit from the fact that 
the respective non-bonding orbitals are mainly centred on 
the metal atom with their centre-of-charge shifted away 
from the bonding region (see FigQ. Of course, an anal¬ 
ogous scheme can be applied to class II molecules. The 
classification of molecules according to [3] is based on 
electronic structure arguments. Class I molecules have 
non-bonding orbitals which appear due to interference 
(hybridization) of atomic valence one-electron wavefunc- 
tions (e.g. s-p or p-d hybridization as it is shown in 
Fig§. Non-bonding orbitals in Class II molecules ap¬ 
pear due to symmetry reasons, like in the HI’*' molecular 
ion [19]. We empasize, however, that even in the case of 
a valence electron undergoing a transition between non¬ 
bonding orbitals, quasidiagonality of the FC matrix can 
only be expected, but not guaranteed. Counteracting are, 
for instance, the Renner-Teller or Jahn-Teller effects, 
which are absent for diatomic molecules, but can consid¬ 
erably influence the equilibrium structure and vibrational 
frequencies of polyatomic molecules. For example, for 
the molecule MgOH the equilibrium structure of the elec¬ 
tronic ground state is linear, whereas the one of the ener¬ 
getically lowest electronically excited state is bent [3D] , so 
that a quasi-diagonal FC matrix can hardly be expected 
in this case. Typically, FC factors are extremely sensitive 
to relative displacements in equilibrium distances of the 
electronic states. Calculation of FC factors with an accu¬ 
racy required for a reliable prediction of laser coolability 
is challenging even for diatomic molecules. Nevertheless, 
even a 10% accuracy in calculations of FC factors, which 
can be conventionally reached in modern quantum chem¬ 
ical approaches for large (bigger than 0.5) FC factors, al¬ 
lows to considerably narrow down the initial selection of 


candidates, whereas simple theoretical concepts primar¬ 
ily provide guidelines for deducing molecular candidates. 

Linear triatomic molecules Examples of linear 
molecules, for which quasidiagonal FC matrices may 
be expected, are calcium monohydroxide (CaOH) and 
calcium monoisocyanide (CaNC). CaOH and CaNC have 
the advantage of featuring a linear equilibrium structure 
in the ground and essentially also in the energetically 
lowest excited doublet state, thus allowing direct adap¬ 
tation of the cooling scheme used for diatomic molecules. 
Besides, CaNC possesses a large electric dipole moment 
(about 6 Debye, see Table I in |3T]), which renders this 
molecule especially attractive for the combined scheme 
of Stark deceleration and subsequent (or simultaneous) 
laser cooling. We calculated molecular parameters and 
FC factors for vibronic transitions between ground and 
excited electronic states. The results are summarised in 
the Table I in |3T], whereas FC factors together with 
vibrational levels are graphically represented on Figj^ 
We considered three largest FC factors as typical cooling 
experiment of diatomic molecules suggest one main 
pump transition and two repumpers as e.g. in |12] . To 
estimate the Doppler temperature for both CaOH and 
CaNC we calculated the natural fluorescence lifetime t 
of the excited states (see Table HI in [21]). Taking r 
being approximately equal to 20 ns, we obtain a Doppler 
temperature Td for both CaOH and CaNC according 
to formula Td = H/(2kBT) where Lb is the Boltzman 
constant and h = /i/(27r) the reduced Planck constant. 
Td is then about 1 mK, which is considerably lower 
than the typical temperature which can be reached in 
the framework of e.g. buffer-gas cooling. 

Computational details: Input files for calculations 
and additional data are also provided in [^. The mol- 
PRO program package [33] was used for electronic struc¬ 
ture calculations. Subsequent computation of FC factors 
in the harmonic approximation were performed with our 
hotFCHT code |33H35] . Basis sets of quadruple-zeta or 
triple-zeta quality augmented with polarisation valence 
basis functions (def2-QZVPP and def2-TZVP) were used 
for all atoms. All electronic structure calculations started 
with the closed-shell Hartree-Fock calculations of the 
ground-state configuration of the singly charged molec¬ 
ular cation to obtain an initial guess for the molecular 
orbitals. Then calculations in the framework of complete 
active space self-consistent field (CASSCF, [3^ ) were 
performed for the neutral molecule, followed by configu¬ 
ration interaction (Cl, [27]) calculations. The FC factors 
for transitions between the ground vibrational state of 
the first electronically excited state and vibrational states 
of the electronic ground state are estimated in all cases. 
A few calculations with different level of accounting for 
electronic correlations and basis sets were used to check 
the stability of the computed FC factors. Details of the 
calculations can be found in [^. We emphasize that in 
all cases the sum of three largest FC factors is always 
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larger than 0.9 for both CaNC and CaOH. This clearly 
indicates the potential of these molecules (and, although 
with exceptions, the general class MOH and MNC, where 
M is an alkaline earth metal Be-Ra) for direct cooling 
with lasers. 

Non-linear polyatomic molecules Examples of non¬ 
linear polyatomic molecules having a valence electron 
on non-bonding orbitals are CaCHa and MgCHa. The 
scheme of arriving at these molecules is analogous to the 
one in the previous paragraph, with the obvious differ¬ 
ence in the replacement of the pseudohalogen by methyl. 
The calculation results are summarized in Table II of 
m- Both molecules have quasi-diagonal FC matrices 
for transitions between electronic ground and first ex¬ 
cited state (see Table II in [21]). The important fea¬ 
ture of these molecules is a quite small Jahn-Teller effect 
(see the data of the experimental measurements of Jahn- 
Teller parameters in |55| for MgCHa), which also favours 
a quasi-diagonal FC matrix. The computational scheme 
was essentially analogous to the one used for CaOH and 
CaNC (same basis sets, calculations of the closed-shell 
cation at first stage followed by MCSCF and Cl calcula¬ 
tions, see ini). As detailed and high-precision electronic 
structure calculations are beyond the scope of the cur¬ 
rent letter, we used a simplified scheme by considering 
only the totally-symmetric modes (Csv molecular sym¬ 
metry group) in estimation of FC factors rather than 
solving the full vibronic coupling problem. To check 
the influence of the non-totally symmetric modes, cal¬ 
culations for MgCHa were performed without symme¬ 
try restrictions and then FC factors were calculated in 
the harmonic approximation. Although the calculated 
molecular parameters are in reasonable agreement with 
experimental measurements (see Table II), the displace¬ 
ments of bond lengths M-C (M=Ca,Mg) predicted from 
our calculations are larger in absolute value than the ex¬ 
perimentally measured ones (see Tables VITX in [21)1. 
Nevertheless, we obtained a sum of FC factors for three 
vibronic transitions with largest FC factors exceeding 0.8 
for MgCHa even when non-totally symmetric modes are 
accounted for. Thus one can hope for even larger values 
of FC factors in experiment. It should be emphasized 
that we primarilly consinder the sum of the three larges 
FC factors, because the value is (according to our ob¬ 
servation) more stable to variations in relative displace¬ 
ments than individual FC factors. In cases where only 
totally-symmetric modes are accounted for, the symme¬ 
try blocks are not perfectly separated due to some resid¬ 
ual symmetry breaking in the excited state structures. 
Nevertheless, we obtained a sum of the three largest FC 
factors exceeding 0.9 for both MgCHa and CaCHa. A 
detailed computational investigation of Jahn-Teller vi¬ 
bronic coupling in these systems shall be left for future 
studies. We note in passing that the quite large elec¬ 
tric dipole moment computed for CaCHa indicates the 
potential for Stark deceleration. Doppler temperatures 


for MgCHa and CaCHa can be estimated analogously to 
linear molecules, using calculated lifetimes of the excited 
states from Table III. Our order-of-magnitude estimate 
of Td is the same as for linear molecules, about 1 mK. 

Chiral molecules Cooling and trapping of chiral 
molecules would mean a crucial improvement for spec¬ 
troscopic measurements of a number of highly interesting 
effects e.g. parity violating energy differences between 
enantiomers due to the fundamental weak interaction 
PSI150] . Using the scheme from the above paragraphs 
it is also possible to identify chiral molecular structures 
having a quasidiagonal FC matrix. One way is isotopic 
substitution to a known molecule |29i m\- The chiral 
methyl group CHDT is well established in chemistry (see 
review [32] for general aspects as well as ref. [33] on parity 
violation in CHDTOH). We calculated FC factors for the 
molecule MgCHDT, using the same method for electronic 
structure calculation as in the case of MgCHa, without 
accounting for molecular symmetry. The resulting FC 
factors we obtained are as follows: (0.64, 0.18, 0.04) with 

equal to 0.86. Thus FC factors are predicted to be only 
slightly smaller than for MgCHa. This modest change is 
excellent news for attempts to obtain (ultra)cold chiral 
molecules for the search of violation of fundamental sym¬ 
metries. 

Conclusions We have outlined a rational approach to 
identify polyatomic molecules that appear to be promis¬ 
ing candidates for direct Doppler cooling with lasers. Ex¬ 
plicit numerical calculations for structures and Franck- 
Condon factors of selected representatives indicate high 
potential for laser-cooling of such molecules to even open 
up the third spatial dimension for ultra-cold molecules 
generated by direct Doppler cooling with lasers. 
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FIG. 1: (Color online) Surfaces of constant value (iso¬ 

surfaces) of the non-bonding one-electron wavefunctions (or¬ 
bitals) of the valence electron for CaCHs. The singly occupied 
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for negative signs of the orbital; the iso value is 0.055 
whereas the SOMO for the excited state configuration is in 
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yellow for negative signs of the orbital; the iso value is 0.040 
It is seen that the electronic centre-of charges are 
shifted outside of the bonding region and that the ground- 
state SOMO is mainly a mixture of s and p wavefunctions 
centered on Ca nucleus (s-p - hybridised orbital), whereas 
the excited state SOMO is rather p-d - hybridised. 
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FIG. 2: (Color online) Schematic graphical representation of the vibronic transitions with largest FC factors for CaNC and 
CaOH molecules. On coordinate axes the normal coordinates Qi and Q 2 are shown for both ground and excited electronic states. 
Qi correspond mainly to Ca-X stretching AR — R — Re, where X=OH, NC, Re is the equilibrium distance Ca-X from Table 
I in m- Coordinate <52 is approximately equal to 9-9e, where 9 is the bending angle and 9e is its equilibrium value according 
to Table I in [^. The approximate displacement vectors for Ca and CN and OH groups are shown next to corresponding 
coordinate axes. All parameters, including the energy E, are in arbitrary units (for the sake of clearer representation). On 
both plots the vibrational quantum numbers and overtones are shown next to the corresponding arrows, e.g !§ denotes the 
vibronic transition between the ground (0) vibrational state of the excited electronic state and the vibrational state with the 
first mode (1) (Ca-X stretching mode) being doubly excited (2) in the ground electronic state. 
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